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ABSTRACT 

In  this  thesis,  the  gain  of  the  magnetic  amplifier  is  derived  theo- 
retically as  a  function  of  different  parameters  of  the  circuit.   An  ex- 
perimental circuit  of  a  series-connected  magnetic  amplifier  with  a 
resistive  load  is  set  up  to  verify  the  theoretical  results. 

Starting  from  Kirchoff's  Voltage  Law  and  Ampere's  Law  applied  to 
the  circuit,  with  the  polynomial  representation  of  the  magnetization 
curve  of  the  core  material,  a  set  of  equations  for  currents  and  fluxes 
is  obtained.   These  equations  contain  nonlinear  terms.  The  Poisson 
perturbation  method  is  used  to  solve  the  simultaneous  nonlinear  equa- 
tions of  fluxes. 

Successive  approximations  of  core  fluxes  and  currents  can  be  made 
to  get  the  theoretical  results  as  close  as  required  to  the  experimental 
results.   Only  the  first  approximation  is  computed  in  this  thesis. 

The  only  difficulty  one  can  encounter  when  using  the  perturbation 
method  is  the  lack  of  knowledge  of  the  region  of  convergence.  The 
smallness  of  the  coefficients  which  are  negative  of  the  nonlinear  fac- 
tors in  this  particular  application  contributes  to  the  rapid  diminution 
of  the  contribution  from  successive  approximations. 
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I-Introduction- 

In  most  analyses  of  magnetic  amplifiers,  the  steady  state  and 
transient  responses  are  obtained  by  piecewise  linear  methods  or  by  the 
assumption  of  a  direct  current  component  and  a  sine  wave  of  carrier  fre- 
quency. 

The  basis  of  operation  of  the  magnetic  amplifier  depends  on  the  non- 
linear characteristics  of  the  core  material.  This  nonlinearity  makes  the 
mathematical  analysis  tedious. 

This  thesis  presents  a  quantitative  analysis  of  the  sinusoidal  re- 
sponse of  a  series  connected  magnetic  amplifier  expressed  in  dimensionless 
form. 

The  method  used  to  solve  the  nonlinear  system  equations  is  I'oisson's 
method.   Kirchhoff's  voltage  law  and  Ampere's  line  integral  law  are  ap- 
plied to  the  circuit  to  set  up  its  fundamental  equations.  The  magnetiza- 
tion curve  of  the  core  material  is  represented  by  a  general  polynomial. 
The  problem  is  then  reduced  to  solving  the  simultaneous  nonlinear  equa- 
tions of  the  fluxes  in  both  cores.   From  thats  successive  approximations 
of  core  fluxes  and  currents  can  be  obtained.   But  only  the  first  approxi- 
mation has  been  done  in  this  thesis. 

The  purpose  of  this  paper  is  to  treat  the  general  problem  of  the 
gain  in  magnetic  amplifiers.   A  general  type  series  connected  magnetic 
amplifier  shown  in  Fig.  1  is  chosen  for  study.   A  biased  sinusoidal  vol- 
tage is  used  as  input;  and  a  power  supply  of  higher  frequencys  as  carrier. 
The  output  is  therefore  a  modulated  signal  subjected  to  harmonic  dis- 
tortions of  both  signal  and  carrier  frequencies.  The  waveform  of  the 
envelope  of  the  modulated  output,  i.e.  the  fundamental  component  and 
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harmonic  contents  of  signal  frequency  has  been  analyzed.  The  following 
results  have  been  obtained  in  dimensionless  form  and  checked  by  experiment: 

a)  General  solutions  for  fluxes  and  currents 

b)  Magnitude  and  phase  angle  of  fundamental  gain  as  a  function  of 
the  system  parameters 

c)  Harmonic  distortions  of  the  response 

Poisson's  method  is  most  valuable  in  cases  where  Poisson's  series 
for  the  flux  converges  rapidly,  and  the  core  material  used  has  a  gradual- 
ly varying  incremental  permeability.  The  problem  of  d.c.  controlled  magne- 
tic amplifier  can  be  analyzed  by  this  method  by  considering  it  as  a 
special  case  of  the  present  general  analysis. 

The  present  analysis  has  the  following  limitations: 
I-If  Poisson's  series  converges  slowly,  the  work  will  be  very  labor- 
ious. 
Il-If  Poisson's  series  diverges,  this  method  will  no  longer  be  valid. 
III-The  core  material  should  have  a  gradually  varying  incremental 

permeability  such  that  it  can  be  suitably  represented  by  a  poly- 
nomial of  few  terms. 


II-Theoretical  Analysis   - 

1.  Assumptions. 

The  magnetic  amplifier  circuit  used  in  this  analysis  is  that  of 
a  series-connected  magnetic  amplifier  shown  in  Fig.  1  with  a  resistance 
load.   The  following  assumptions  are  made: 

a)  Both  cores  have  the  same  electrical  and  magnetic  prop- 
erties, 

b)  Effects  due  to  eddy-current  and  hysteresis  losses  in 
the  core  are  neglected, 

c)  The  magnetization  characteristics  can  be  suitably  re- 
presented by  a  polynomial 

2.  Fundamental  equations. 

Applying  Kirchhoff's  voltage  law  to  the  circuit  in  Ftg.  2 


tjhcs.    t   JV,  tlj3f>  (i) 

t'C  e  it 


'A 

Applying  Ampere's  Law  to  the  magnetic  circuits  of  the  cores  (core  a 
and  core  b)  (Fig.  2). 

Fa  -  Wx  lx  +  N2  i2  (3) 

Fb  -  Nl  ll  "  N2  l2  <4> 

The  magnetization  characteristics  of  the  magnetic  saturable  reactors 
is  represented  by  a  polynomial  containing  a  sufficient  number  of  terms 
to  secure  the  exactness  of  fit.   In  general: 


F  -  !L   k  ^  (5) 

a   ~—  n  a  ' 


where  n's  are  odd  integers  only,  due  to  the  skew- symmetry  of  the  magnetiza- 
tion characteristics. 

The  6  above  equations  constitute  the  fundamental  system  equations  of 
the  magnetic  amplifier  circuit  shown  in  Fig.  1.  These  will  be  solved 
to  determine  the  transient  andsteady  state  of  currents  i  and  i_. 

Substituting  aquations  (5)  and  (6)  into  equations  (3)  and  (4)  re- 
spectively, then  solving  for  i  and  i_,  in  terms  of  'p 


u  =  -    i  «  f  *£  k  <T  (7) 

i,  -  —  I  it  0*  -  £  *    K\  (8) 


Combining  corresponding  terms  of  the  polynomials  in  the  above  equa- 
tions, 

i       ^    * 

Substituting  equation(9)  into  equation  (1)  and  equation  (10)  into 
equation  (2), 

E?,   ,.   d  :   ■^■",:-.  >  t   ;;".   '-''    < ;  (12) 

Let 

0a  +  0b  -  20x  (13) 

0a   -  0b  =  20y  (14) 

we  may  then  write, 

(15) 
£'*,.  ..  cj  (16) 
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U 


7<&JUUUlJ 


< 
CD 
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t 0$TTPs  r^WTO^|_ 
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Q. 
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where  the  j8s  are  odd  integers  greater  than  unity. 

Substituting  equations  (13),  (14),  (15)  and  (16)  into  equations  (11) 

and  (12),  then  rearranging  the  results,  the  following  nonlinear,   simul- 


taneous equations  for  0  and  0  are  obtained: 

x     y 

ctfc 


volt/turn   (17) 


-«  '-  , 

=  l- 

T   Vr 

)      volt/turn 

• 

(18) 

here 

T  = 

N 

second 

(19) 

i  ft' 

second 

(20) 

V.   , 

' 

volt /turn 

(21) 

V 

volt/turn 

(22) 

v 

'■ 

volt /turn 

(23) 

I  = 

dimensionless 

(24) 

1 

p 

dimensionless 

(25) 

Rearranging  the  terms, 


1      I  ■* 


(26) 


f.  •'-.  4, 

,*  .       ... 


A 


(27) 


where  -*  is  a  dimensional  factor  which  is  introduced  in  these  equations 

in  order  to  make  the  coefficients     and     dimensionless.   Numerically 

2 
<>    =  1  and  has  dimensions  of  turns   ,  ohm.   By  examining  equations 

(13)  and  (14)  and  the  circuit  diagram  in  Fig.  2,  it  is  found  that  0  is 

the  flux  principally  produced  by  the  control  current  and  0  is  the  flux 


produced  by  the  output  current. 
In  equation  (18): 

volt/turn    (18) 
The  left  hand  side  term  is  the  generated  counter  emf  along  the  output 
circuit,  and  the  first  term  of  the  righthand  side  is  the  linear  part  of 
the  voltage  drop  of  the  fundamental  component  of  the  carrier  frequency; 
the  second  term  is   the  applied  carrier  voltage  and  the  last  term  is 
the  voltage  drops  of  all  harmonics. 

Equation  (17)  has  a  corresponding  physical  meaning  for  the  input 
circuit. 

From  equations  (13)  and  (14), 


0     =  0     ■?•  0  weber  (28) 

a    x    y 

0,    ■  0     -  0  weber  (29) 

b    x    y 

Substituting  equations  (29)  and  (30)  into  equations  (9)  and  (10), 

the  solutions  of  the  current  responses  in  terms  of  0     and  0'  are  obtained: 

x      y 

lit.-., 
I     ■-        f.  .- 


i,  =  -  f  (30) 


I  -  i  '  (31) 

! 

If  the  order  of  the  polynomial  in  the  above  equations  is  allowed  to 
increase  indefinitely,  equations  (30)  and  (31)  become  infinite  series. 
The  convergence  of  this  infinite  series  is  proved  in  Appendix  A. 

Equations  (30)  and  (31)  for  the  current  responses  can  be  utilized 

only  if  simultaneous  nonlinear  equations  (17)  and  (18)  can  be  solved  for 

0     and  0   .   In  these  nonlinear  equations,     ,  the  coefficient  of  the 
x      y  *     \\     * 

nonlinear  factor   f  ,  and     ,  the  coefficient  of  the  nonlinear  factor 
f^  are  much  less  than  unity  by  virtue  of  the  low  internal  resistances  of 


FIG.  2-   SERlES-CONNECTEb     MAGNETIC  AMPLIFIER 
(     Theorefical  Circuit     ) 


both  control  and  gate  circuits.  In  case  that  the  absolute  values  of  the 
nonlinear  terms  and  are  less  than 

that  of  the  linear  terms  (-  1 0  t  \,       A**  .*it  +  V  '^  and  (  -  ^*i~c) 

respectively,  equations  (17)  and  (18)  can  be  assumed  to  be  quasi- linear. 
Since  the  nonlinear  terms  appear  with  small  coefficients  :\  and  , 
the  solutions  of  equations  (17)  and  (18)  can  be  expressed  in  the  form  of 
Poisson's  series: 

0X  -   ^       4   ,  1  (32) 

0y  -    :  ^  •:.  ,.  (33) 

where  v    and  :   are  given  by  equations  (24)  and  (25)  respectively. 
Rearranging  equations  (32)  and  (33), 


|    '  V:.       -,    .    j         : 


*y   -  if  -    ,  v  (35) 

Substituting  equations  (32)  and  (35)  into  equations  (17),  (33)  and  (34) 

into  equation  (18)  then  equating  terms  containing  the  same  powers  of 

and   \   respectively,  the  following  series  of  linear  differential  equ?  - 

tions  are  obtained: 


(36) 
(37) 
(38) 
(39) 
(40) 
(41) 


i  ■ 

i  4.    *  Id 

&  A           .:    d 

" 

■ 

'■','.                       C; ' , 

i   4,  <; 

'    .. 

*\    .,'•.       •     v 

v\   '••  • 


t 


'  • 


In  the  above  equations,  f„   and  f_    are  the  partial  derivatives  of 

f  with  respect  to  0  and  0  . 

x      y 

Equations  (36)  and  (41)  are  linear.  The  first  set  of  equations  (36) 

and  (37)  can  be  easily  solved  for  0   and  0     .  This  solution  is  substituted 

xo      yo 

into  equations  (38)  and  (39)  which  can  then  be  solved  for  0  ,   and  0  .. 

7  xl      yl 

This  solution  is  substituted  into  equations  (40)  and  (41)  to  solve  for 
0  ?  and  0  «,  and  so  on.... In  general  the  solution  of  kth  set 
leads  to  the  solution  of  the  (k+l)th  set. 

In  all  cases  in  which  convergence  of  Poisson's  series  given  by  equa- 
tions (32)  and  (33)  exists,  the  exact  solution  can  be  approached  by  reitera- 
tion of  successive  substitutions.   The  convergence  of  Poisson's  series  is 
proved  in  Appendix  B.  The  additional  term  obtained  by  a  further  sub- 
stitution is  smaller  than  the  previous  one  and  ultimately  further  sub- 
stitution leads  only  to  a  negligible  correction.   In  many  practical  cases 
the  respective  first  term  0   and  0   of  equations  (32)  and  (33)  alone 
are  sufficiently  close  to  the  exact  answer.   Such  a  solution  is  the  first 
approximation  solution.   In  the  first  approximation, 

0   =0  (42) 

X     xo  ' 

0   =0  (43) 

y      yo 

A  closer  solution  can  be  obtained  by  taking  the  first  two  terms  of 
the  series  (32)  and  (33),  and  the  solution  is  called  the  second  approxi- 
mation.  In  the  second  approximation, 

0   =0   +  0  .  (44) 

x     xo    xl 

0   =0   +  0  .  (45) 

y     yo    yl 
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If  a  more  exact  solution  is  wanted,  more  terms  in  the  series  (32) 
and  (33)  should  be  taken,  and  the  solution  is  denoted  as  kth  approxi- 
mation, in  which, 

*«-  4u>  +  fA,  +    -   •  -  ,       <*6> 

t     -  <k       .         fa  (47) 

3.   First  approximation  theory. 

For  a  practical  application  of  the  theory  derived  in  the  previous 
sections,  and  to  simplify  the  equations  without  affecting  seriously  the 
exactness  of  the  polynomial  to  the  F,  0  curve,  it  can  be  assumed  that  the 
magnetization  characteristics  of  the  core  material  can  be  approximated 
by  a  three  term  polynomial  within  its  operating  range, 

3       5 
F  »  k  0  +  k,  0  +  k_  0       ampere-turn  (48) 

and  that     and    are  so  small  such  that  equations  (32)  and  (33)  can 
be  approximated  by 


0*0  weber     (49) 

x    xo 

0  ■  0  weber     (50) 

y    yo 

Evaluating  0   and  0   from  equations  (36)  and  (37) „  and  substituting 
them  into  equations  (49)  and  (50),  the  solution  of  fluxes  of  first  approxi- 
mation is  obtained, 

0  =  0  .  +  0,  sin  (.   t  -  0.)  weber     (51) 

x    ob     Lm  1 

0  m   0.   sin(  UJ  t  -  9_)  weber     (52) 

y     Zm  e. 


where 


«,b  ■  sji,  »eber      <53) 


11 


N  E 

0,      -    -, 7£— •— >^-"1-1  "/  ■ — "-  <2  weber     (54) 

lm    kl   Rl  V  {l   +  (    ^ 


0 


2 2m  weber     (55) 

k  2R2  V     (1  +  RT  )2  +  (  S  )2' 
2R 


2 

Gj  -  tan  degree    (56) 

0~  ■  tan  degree    (57) 

k   R 
J\  =  radian/sec   (58) 

Nl 

kl  R2 
U* .  = radian/sec   (59) 

A 

Substituting  equations  (51)  and  (52)  into  equations  (30)  and  (31)and 
letting  n's  be  the  odd  numbers  up  to  5,  the  current  responses  of  the  first 
approximation  attain  the  following  terms: 


l.  -   I   +1   sin   (  „_t   -  9.  )  +  I  ,  cos  2(  t  -  0.  ) 

1    3   -  •  8       ■  ' 

+  I  sin  3(  It  -  e  )   +  I  cos  4(  t  -  e  )         (60) 

f 

+  I  sin  5(  ,t  -  6  )   +  %  sin  ( rt  -  0  )  +  I 

+  1  cos  2(  .  t  -  6  )   +1   sin  3(  .t  -  0  )  |  cos  2(uJt  -0  ) 


+  \l     +  I  sin(<  t  -  0  )    cos  4(^t  -  0  )  ampere 


i0  -  !l  +  I  sin(/  t  -  0  )  +1  cos  2(ut  -  0  )  +  I  sin  3(.ftt  -  0  ) 

+  I  cos  4(^t  -  0)|  sin  (xt  -  0  )   +  I  I  +  I.  sin(^t  -  ©) 


+  I  cos  2(.,t  -  0  )\   sin  3(xt  -  0  )  +  I.  sin  5(..t  -  ©  )  ampere 

(61) 
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where  the  I*s  are  given  in  Appendix  C.   If  we  consider  the  envelope  of 
the  modulated  signal  of  only  the  fundamental  carrier  frequency  across 
the  load  resistance,  the  output  voltage  will  be: 

eQ  =  R   I  sin  (  <  t  -  0  )  +  I  cob   2(  <t  -  0) 

+1  sin  3(.'.it-e,)  +  I  cos  4(^-0)  (62) 

The  fundamental  voltage  gain  is  then: 


v      E. 
lm 

dimensionless  (63 

and  the  percentage  harmonic  distortion: 

%   2nd  Harmonic  distortion  » 

hi 

-  X   100 

(64) 

%   3rd  Harmonic  distortion  ■ 

T13 
Xll 

-  x  100  , 

etc (65) 

The  phase  difference  between  input  voltage  and  the  fundamental  component 
of  the  output  voltage  is  evidently  0  . 
4.   Fundamental  voltage  gain. 

From  the  value  of  L.  defined  in  equation  (155)  if: 

E,   -  .   'tl~  /  3    I  kj  \         volts     (66) 

eio  -  ~^~  JJT^\      - volts    (67) 

e20  *    *Sl*l  JTTKT  volts        (68) 

A  -  1,0%    ,3    ZL     ~     dimensionless   (69) 
the  fundamental  voltage  gain  derived  in  equation  (63)  can  be  expressed 


as: 


p  v*- 


K *   l-(Ei  -h±_  .  lil!        (70) 
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where  «r», ..   and  -uitJ   were  defined  in  equations   (58)  and  (59). 

The  equation  (70)  shows  that  this  gain  is  a  function  of  six  dimension- 
less  parameters; 

Eb    El   E2   \       ^    and   -rL 
Ebo  E10   E20   7*2^  -"■* 

5.  Amplitude  of  voltage  gain  as  a  function  of  the  dimensionless  para- 
meters: 

(a)  Gain  as  a  function  of  bias  - 

If  everything  is  considered  constant  except  E./E,  ,  the  gain 

b  bo 


may  be  expressed  as 
K  - 


febo     L 


where 


.  /E-x  \p\     |.l^\  I'l^Vii"  f  (-"  1    (i*!i?  \ 


and 

,v"  <■     ■    IjT-  :  nu,'       ST 

The  theoretical  characteristics  gain  curve  as  a  function  of  bias 
is  plotted  in  Fig.  3. 

When  the  bias  is  zero,  the  signal  voltage  generates  a  corresponding 
flux  in  one  direction  in  one  half  cycle,  then  in  the  other  direction  the 
following  half  cycle.   The  envelope  of  the  modulated  output  is  dominated 
by  even  harmonics.   The  fundamental  gain  is  zero.   Increasing  the  bias, 
which  prevents  the  flux  generated  by  the  signal  voltage  from  being  of  the 
reverse  direction,  results  in  an  increase  of  fundamental  gain  and  decrease 
of  even  harmonic  distortion.   If  the  bias  is  further  increased,  the  second 
order  term  in  the  bracket  of  equation  (71)  becomes  significant  and  K  will 
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FIG  .  3    -  GAIN        V/  S 

(      Theoretical    ) 
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BIAS 


decrease.   Therefore  the  curve  of  the  gain  will  bend  down.   Physically, 
it  means  that  the  bias  has  caused  saturation. 

(b)  Gain  as  a  function  of  carrier  voltage- 

The  amplitude  of  the  carrier  voltage  will  produce  the  same 
effect  as  the  bias.   The  characteristic  gain  curve  is  similar  to  that  of 
Fig.  3. 

(c)  Gain  as  a  function  of  input  voltage  - 

The  gain  is  independent  of  the  input  voltage  as  long  as  this  input 

voltage  is  small.   When  the  input  voltage  increases,  the  third  term  In 

2   2 
the  bracket  of  equation  (70)   E./Eir.   becomes  significant  in  comparison 

with  the  summation  of  the  other  terms  in  the  bracket,  therefore  the  gain 
will  decrease.   Physically  the  flux  swing  caused  by  a  large  signal  vol- 
tage covers  a  large  portion  on  the  nonlinear  magnetization  curve.   More 
harmonics  are  thus  generated  while  the  fundamental  component  does  not 
increase  proportionally  with  the  input. 

If  the  dimensionless  input  voltage  E./E.    is  considered  as  the  only 
variable  in  equation  (70),  the  gain  may  be  expressed  as: 

e,  •* 


(74) 


where 


and 


\* 


(75) 


,  (76) 


Equation  (74)  is  plotted  in  Fig.  4. 

(d)  Gain  as  a  function  of  signal  frequency  - 

If  the  dimensionless  signal  frequency   /      is  considered  as 
the  only  variable  in  equation  (70),  the  gain  will  be  expressed  as: 

16 


FIG. 4-        GAIN      V/  S 
(    Theoretical    } 
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IA 


INPUT 


K 


where 


I   - 


A 


(77) 


«* 


^&Mi)H  »-^Uc-Sj^T} 


JV. 


and 


i  +  St  V4, 


cl/]l/i 


(78) 


M" 


j    P.    \* 


i  i* 


2s  Y 

aft* 


UBurT 


(79) 


(e)  Gain  as  a  function  of  carrier  frequency  - 

The  carrier  frequency  has  a  similar  effect  upon  the  gain. 
Plotted,  the  gain  characteristic  curve  will  be  similar  to  the  curve  in 
Fig.  5. 

(f)  Gain  as  a  function  of  load  resistance  - 

If  the  dimensionless  load  resistance   rt/2r  *s  tne  only 
variable  in  equation  (70),  the  gain  will  be: 

«^h.\  r 


K 


'L 


(80) 


«..  \* 


where 


(81) 


and 


h 


( 


<-ifcr-i£n>^'] 


-i 


(82) 


For  a  definite  value  of  w/„  ,  the  gain  curve  is  plotted  in  Fig.  6. 
According  to  the  gain  equation  (Eq.  80),  the  gain  would  reach  asymptoti- 
cally a  finite  value  as  load  resistance  increases  to  infinity.   However 
physically  the  gain  is  expected  to  approach  zero  as  load  goes  to  infinity. 
When   R  goes  to  infinity,  the  gate  circuit  acts  as  an  open  circuits, 
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GAIN       V/S        SIGNAL       FREQUENCY 
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and   1      €    AC*~   U/t  s     tl 


since  t     is  independent  of      which  drives  the  core  into  saturation. 

6.  Phase  angle  - 

The  phase  difference  8-  between  input  voltage  and  the  fundamental 
component  of  output  appears  implicitly  in  equation  (62).   Let  L.ft  be 
the  linearized  inductance  of  the  input  winding, 

L1Q  «  N*  /  k.  henry  (83) 

From  equations  (56)  and  (58) 

0.  »  tan    JO  degree  (84) 

Ri 

The  phase  difference  between  input  current  and  output  current  is 
negligible. 

7.  Harmonic  distortion. 

(a)  In  the  output  voltage  equation  (Eq.  62),  the  first  term  R  I 
represents  the  fundamental  component,  the  second  term  R.I-.^  represents 
the  second  harmonic  content,  the  third  term  R.I,-  represents  the  third 
harmonic  content  etc. 

Only  the  second  harmonic  distortion  is  considered  in  this  presenta- 
tion. 

Let  the  second  harmonic  generation  be  defind  as; 

K2nd  -  Il2  RL  (85) 

or  as  a  percentage  of  the  fundamental  gain: 

2nd 

Z  2nd  harmonic  distortion  =  — (86) 

K 

K  is  the  fundamental  gain  given  by  equation  (70).   If  I.-  is  substituted 

O  A 

by  its  value  given  in  equation  (156)  into  equation  (85),  K    may  be  ex- 
pressed as: 

20 


2 


FIG.- 6    - 


GAIN        V/S       LOAD      RESISTANCE 

(    Theoretical  ) 
21 


K2nd, ■•■       I,. a  a      ~\     -    -  (87) 

where  R*        "°  ^ 


2nd 


1       i 

;    *  (88) 


' 


(b)  Effect  of  bias  on  K  and  K2nd  - 

From  equations  (70)  and  (87)  and  from  Fig.  3  and  7,  it  is  seen 
that  as  bias  approaches  zero,  K    approaches  a  finite  value  and  K  approach- 
es zero.   In  the  limit,  the  absence  of  the  fundamental  frequency  and  a 
maximum  second  harmonic  in  the  output  are  expected.   In  fact,  when  the  bias 
is  zero,  odd  harmonics  are  absent. 

(c)  Second  harmonic  generation  as  a  function  of  bias  - 

If  the  dimensionless  E,  /E,   is  considered  as  the  only  variable 

D        DO 

equation   (87)   becomes: 

K2nd     -      .<2nd     I     ,..     ?ni         h     V   1     •  (89) 

where  -i 

.,  awf!i\(|     ;    t    ■■      :;--..,.:  .  ■  i 


,.      £   r]    I    ...    '  (90) 


and 


»»  .      -     e.  .  <9i) 


i 


j 


The  second  harmonic  generation  as  a  function  of  bias  is  plotted 
in  Fig.  7. 
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SECOND     HARMONIC      GENERATION 


PLOTTED     AGAINST      BIAS    (Theoretical) 
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Ill-Experimental  Verification  - 
1.   Experimental  setup  - 

(a)  Experimental  circuit  - 

The  circuit  used  for  this  experiment  is  illustrated  in 
Figs.  1  and  2.   Two  identical  saturable  magnetic  reactors  are  connected 
in  series.  An  alternating  current  voltage  source  biased  by  a  direct 
current  voltage  source  was  used  to  provide  the  input  signal  to  the  cir- 
cuit. An  a  c  voltage  source  of  higher  frequency  was  used  as  carrier 
supply. 

(b)  Equipment  and  instrumentation  - 

The  cores  used  were  a  matched  pair  of  core  #50106  (HYMU80, 
Magnetic  INC.  Products).  An  Hewlett-Packard  audio  generator  type  HP  202A 
was  used  as  signal  generator.   Between  the  signal  generator  and  the  in- 
put circuit  a  cathode  follower  was  connected  so  that  negligible  impedance 
has  bean  added  to  the  control  circuit.   The  cathode  follower  used  in  this 
experiment  is  illustrated  in  Fig.  18.   The  carrier  supply  for*  the  magnetic 
amplifier  was  taken  from  the  400  cycle  per  second,  120  volt,  2-phase  labor- 
atory bus,  through  a  powerstat  to  control  the  carrier  voltage.  The  bias 
voltage  was  supplied  by  a  6  volt  d  c  source  through  slide  wire  resistor 
to  provide  bias  voltage  desired.  A  slide  wire  resistor  was  connected  in 
series  with  the  bias  source  to  keep  the  resistance  of  the  control  circuit 
constant  for  different  bias  voltages.  The  load  resistance  used  was  an- 
other slide  wire  resistor. 

The  cores  were  wound  with  75  turns  of  AWG#38  wire  for  the  control 
circuit  and  500  turns  of  AWG#41  wire  for  the  gate  circuit.  The  load 
voltage  was  displayed  on  a  Tektronix  Oscilloscope  Type  545;,  using  type  K 
plug-in  amplifier  for  single  trace  and  type  CA  plug-in  unit  for  dual  trace. 
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An  General  Radio  Wave  Analyzer  type  736-A  was  used  to  measure  the  magni- 
tude of  the  second  harmonic  in  the  output  current. 

(c)  Magnetization  curve  of  the  magnetic  reactor  - 

The  magnetization  characteristics  of  the  saturable  magne- 
tic reactor  were  measured  by  using  the  circuit  arranged  as  shown  in  Fig. 8, 
The  display  of  the  X  axis  of  the  oscilloscope  measures  the  magnetomotive 

force  F  and 

x 
F  ■  N^  ampere-turns   (92) 

The  Y  axis  measures  the  total  flux  when  the  values  of  r_  and  C-  are  so 

designed  that 

r  vv     1  ohms  (93) 

"     U,C2 

where  ui  is  the  frequency  of  the  power  supply. 

Since  the  voltage  equation  of  the  secondary  winding  of  the  reactor 

(Fig.  8)  is 

N2"f    -  V2  V-fv'-  (94) 

C2 
According  to  the  relationship  (93),  the  reactance  drop  in  equation  (94) 

can  be  neglected  without  introducing  serious  error,  thus 

The  voltage  applied  to  the  Y  axis  of  the  oscilloscope  is 

V  =  pi—  I     i  dt  (96) 

y     2  J 

Substituting  the  value  of  i~   in  equation  (95)  into  equation  (96), 

r  C 
7     7       V 

0  m       „  —  y  weber         (97) 

N2 

Equation  (97)  shows  that  the  total  flux  is  proportional  to  the  display 
of  the  Y  axis  of  the  oscilloscope.  Thus  the  magnetization  curve  of  the 
saturable  reactor  can  be  obtained  simply  by  calibrating  the  oscilloscope., 
The  magnetization  characteristics  of  the  reactors  used  in  this  experiment 
are  shown  in  Fig.  9.  25 
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Circuit      for       measuring    the     magnetization 
characteristics       of        the 

Magnetic     Reactors 
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(d)  Computation 

The  magnetization  curve  shown  in  Fig.  9  may  be  represented 
by  a  three-term  polynomial  given  in  equation  (48).   Using  numerical 
analysis  method  to  determine  the  value  of  k  ,  k?,  and  k  ,   there  results 

k  =  .177  x  10  amp- turns /weber   (98) 

12 
k2  =  -11.6  x  10  amp-turns/weber   (99) 

18 
k  *  348  x  10  amp-turns/weber   (100) 

The  magnetic  amplifier  has  the  following  constants: 

Nx  «   75      turns  (101) 

N2  =  500      turns  (102) 

R.  =    6.05   ohms  (103) 

R2  «   113      ohms  (104) 
From  equations  (58),  (59),  (66),  (67),  (68)  and  (69),  the  following  con- 
stants are  evaluated: 

■D-e      m     196     radians/second  (105) 

WQ  s   80     radians/second  (106) 

E.   *    2.85  volts  (107) 

DO 

E.   =    2.32  volts  (108) 

lo 

E20  =    6'5   V°ltS  (109) 

A   =  12.8     dimensionless  (110) 

If  the  following  specific  numerical  values  are  used  for  the  computations, 

E2  -  95       volts  (111) 

UJ   »  2«  x  400   radians/second  (112) 

E   -  1.00     volt  (113) 

„fv-   s  2rt  x  20   radians /second  (114) 
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E20 

U/. 

=  31.4 

El 

=  0.43 

Eio 

=  0.64 

Eb 
Ebo 

*  (0.512) 

RL 

=  22 

2R2 

Eb  =  1.46          volts  (115) 

RL  «  5000          ohms  (116) 

the  six  dtmensionless  parameters  which  are  present  in  the  fundamental 
voltage  gain  equation  (Eq.  70)  will  be: 

=  14.4  (117) 


(118) 
(U9) 

(120) 
(121) 

(122) 


and  the  coefficients   P   and  0      of  the  nonlinear  factors  given  by 
equations  (24)  and  (25)  are: 

(\   =  -0.54  x  10"3  (123) 

l(      m  -   5.20  x  10"3  (124) 

From  equation  (70),  the  gain  of  first  approximation  can  be  calculated: 

K  «  24.14  (125) 

(e)  Experimental  results  - 

Fig.  11  shows  the  waveforms  of  the  envelope  of  the  output 
currents  for  different  biases  with  other  parameters  being: 

Ej  -  1.0  volt  (126) 

E2  -  95  volts  (127) 

il  a  2  x  20  radians/second  (128) 

UJ  s  2  x  400  radians/second  (129) 

R  «  5000  ohms  (130) 

Li 
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FIG.  10-       EXPERIMENTAL  SET-UP 


and  the  calibration  of  the  oscilloscope  being  8  volts  per  small  division 
of  the  screen  of  the  oscilloscope. 

With  Efe  »  1.46  volt  for  example,  the  peak  to  peak  of  the  modulating 
fundamental  covers  8.1  small  divisions.  Thus  the  gain  is 

K  =    8A   X  8, -  23.0  (131) 

2  x   V*  x  1.0 

From  equations  (72)  and  (73), 

=  63  (132) 

>j   =  1.48  (133) 

From  equations  (131)  and  (132), 

_K »  0.365  (134) 

Ai 

Since  the  constant  E,   =  2.85  volts  (Eq.  107)  and  the  bias  voltage  E,  = 

1.46  volt, 

Eb  =  1>46  -  0.512  (135) 

E.     2.85 
bo 

By  equations (134) and  (135),  this  particular  experimental  result  can  be 
plotted  with  K/,    as  ordinate  and  E,  /E,   as  abcissa  for  the  experiment- 
al curve  of  gain  versus  bias.  (Fig.  12). 

Fig.  11  (d)  shows  negligible  second  harmonic  distortion,  however, 
in  the  case  of  zero  bias  (Fig.  11  a),  the  modulating  wave  is  dominated 
by  second  harmonic,  which  covers  approximately  1.6  small  division  of  the 
screen  of  the  oscilloscope.  Then 

1.6  x  8     -  4.57  (136) 


(137) 
(138) 
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FIG.  II-     ENVELOPES      OF     OUTPUT         CURRENTS 
FOR  DIFFERENT         BIASES 
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FIG.    12  -    GAIN         V/S         BlAs 
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Thus  for  E,  /EL  =  0 
b  bo 

K2nd    -  4^57   «  0.75  (139) 

<*«**     6.1 

b 

Different  values  of  K  n  /  .,;  n   with  corresponding  values  of  E^/E, 

b  bo 

are  plotted  in  Fig.  12. 

Fig.  13  shows  the  envelopes  of  the  output  currents  for  different 
signal  frequencieswith  other  parameters  being 

E1  -  1.0  volts  (140) 

E2  =  95  volts  (141) 

^  ■  2it  x  400  radians/second  (142) 

Eb  «  1.46      volt  (143) 

RT  =  5000   ohms  (144) 

and  calibration  being   6  volts  per  small  division  of  the  screen  of  the 
oscilloscope. 
From  equation  (78)  and  (79) 

<*a  =  32  (145) 

hi    =0.3  (146) 

Different  values  of  K/x ,-    with  corresponding  values  of  -1/.   are 
plotted  in  Fig.  14. 

Fig.  15  shows  the  envelopes  of  output  currents  for  different  load 
resistances.   An  experimental  curve  is  obtained  in  Fig.  16  with  different 
values  of  K/^c,    plotted  against   IL/2R-. 

(f)  Comparison  of  experimental  results  with  computation 

In  Fig.  11  (Gain  versus  bias),  the  theoretical  curves 
approach  zero  much  more  rapidly  than  the  experimental  ones  if  the  bias 
is  brought  up  into  the  region  of  over-saturation.   A  deviation  is  expect- 
ed in  this  region  where  the  three-term  polynomial  given  in  equation  (48) 
does  not  closely  represent  the  magnetization  curve. 
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(a  )       ^  =    21T    K20 


(  b  )      -O-    =    2  tT    x  30 


(c)      ./*     =    2  TT    x   50 


(d)      A    =  21T    X80 


(e)     H     =    2TT    XI20 


riG.    13-      ENVELOPES        OF       OUTPUT       CURRENTS 


FOR         DIFFERENT  SIGNAL       FREQUENCIES 
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FIG.  14-        GAIN      V/S         SIGNAL         FREQUENCY 
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Similarly  for  the  case  of  the  varying  load  resistance   (Fig.  16). 
For  a  small  load  resistance,  both  equation  (80)  and  experimental  results 
show  a  linear  relationship  between  the  gain  and  the  load  resistance.   If 
the  load  resistance  is  increased,  this  relationship  is  no  longer  linear, 
but  there  still  exists  an  agreement  between  the  theoretical  and  experi- 
mental results.   If  the  load  resistance  is  further  increased  beyond  a  cer- 
tain limit,  agreement  is  no  longer  expected. 

It  becomes  evident  that  if  the  load  resistance  is  increased  beyond 
a  certain  limit,  the  first  approximation  is  no  longer  sufficient  and 
second  or  even  higher  approximation  becomes  necessary.   It  is,  there- 
fore, necessary  to  establish  a  limit  for  load  resistance  beyond  which 
the  first  approximation  is  insufficient. 

It  is  mentioned  in  Appendix  B  that  if  the  first  approximation  0 
is  to  be  sufficient  and  if  the  term  in  equations  (32)  and  (33)  monotonical- 
ly  decrease,  the  additional  term  0  -  should  be  such  that 

|fa|l^iU<     l<V"l  (143) 


or 


■ 


"    "  (144) 

Since   from  equation   (25) 


I3     7~2 

2  N2 


!2 


R2   (   1  *2rT  I  (145) 


the   inequality   (144)  may  be  written  as 

J2  -    1  !   R. 


0  ?m2 

yo  : 


*L    «  IT-1-      -4-       "    l  2R2  (146) 


0-l  •  R2 


This  establishes  an  upper  limit  for  load  resistance  for  which  the  first 
approximation  is  sufficient. 
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(d  )     RL        =  15000   ohms 


FIG.  15-      ENVELOPES  OF        OUTPUT        CURRENTS 

FOR         DIFFERENT         LOAO        RESISTANCES 
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FIG.J6-      GAIN       V/S      LOAD       RESISTANCE 
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FIG.  17-     SECOND    HARMONIC    GENERATION 

PLOTTED    AGAINST     BIAS 


IV-Conclusion  - 

The  present  analysis,  principally  concerned  with  the  treatment  of 
a  series-connected  magnetic  amplifier  including  the  effects  of  a  re- 
sistive load,  is  mainly  based  on  Polsson's  perturbation  method.   The 
lack  of  knowledge  of  the  region  of  convergence  is  the  main  disadvantage 
of  the  use  of  this  method.  This  type  of  difficulty  is  inherent  in  most 
perturbation  methods.   In  this  particular  application  the  coefficient 
of  nonlinear  factor  is  usually  small  and  negative.   The  smallness  con- 
tributes to  a  rapid  diminution  of  the  contribution  from  successive  approxi- 
mations.  In  addition,  its  negative  sign  results  in  a  series  whose  terms 
are  alternately  positive  and  negative.   If  the  absolute  values  of  the 
terms  form  a  monotonic  null  sequence,  the  series  converges.   Further- 
more, the  error  that  results  from  approximating  the  infinite  series  by 
a  finite  number  of  terms  does  not  exceed  the  absolute  value  of  the  first 
term  omitted. 
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FIG. 18 


CATHODE      FOLLOWER 


used  in  the  input  circuit  of  th&    experimental    setup 
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APPENDIX  A 

Proof  of  convergence  of  Equations  (30)  and  (31). 

Since  0  and  0  are  finite,  there  may  be  assigned  a  finite  value 

0  /2,  by  which  both  0  and  0  are  bounded,  i.e. 
o  x      y  * 

0X  4  0Q/2  (147) 

0„  ,  0/2     if   0  >  0  (148) 

y     o  x     y  v   ' 

Substituting  (147)  and  (148)  into  equation  (31),  the  latter  is  reduced 
to 


1 


1,3,5...  l,3,5...r 


i,  ..  rr~  \        k     0  \      n 

2   N0  n     o  f 


2 


^  m 


C)  <«»> 


bmce 


n  '  -  2n  (150) 


i  m  . 

then 

i     i,-J,j... 

i,  <  tt-  k  0n  (151) 

2    N2  no  ' 

n 

Because  the  right  hand  side  of  equation  (151)  is  a  series  known  as 
convergent,  therefore,  the  current  response  equations  for  i-  Is  con- 
vergent too.  The  convergence  of  i  which  is  given  by  equation  (30) 
can  be  proved  the  same  way. 


43 


APPENDIX  B 
Sufficient  conditions  for  the  convergence  of  the  Series  given 
by  equations  (32)  and  (33) 

The  assumption  that  0  and  0  may  be  represented  by  equations  (32) 

*    y 

and  (33)  is  valid  only  if  such  series  do  converge.   Otherwise,  0  and  0 

x      y 

cannot  be  represented  by  these  equations  and  the  present  analysis  is  not 
applicable. 

A  sufficient  condition  for  the  convergence  of  these  series  can  be 
established  due  to  the  fact  that  the  small  parameters     and      are 
negative.   Thus,  odd  power  term  of  equations  (32)  and  (33)  are  of  nega- 
tive sign  and  even  power  term,  of  positive  sign.   There  results  a  series 
whose  terms  are  alternatively  positive  and  negative.   The  sufficient  con- 
dition for  the  convergence  of  such  series  according  to  Leibnitz  theorem 
is  that  the  absolute  values  of  the  terms  form  a  monotonic  null  sequence, 
namely,  for  the  present  case, 


\  ^x(n-l)l  (152) 

;\l  J  10 

\    x  n  i 


and 


l« 


-0v(n-l)  (153) 


y  n 


Furthermore,  the  error  that  results  from  approximating  the  infinite 
series  by  a  finite  number  of  terms  does  not  exceed  the  absolute  value  of 
the  first  term  omitted. 


44 


APPENDIX  C 
Expressions  of  currents   appearing   in  Equations    (80)   and   (61) 


I10-     h.     *2m[l  +  3     ^3      (0^     +I     **l     *\$\j 
N2  kx 


(154) 


c  kc  /rf  ''  -3  rt4          1  M4          _   .2  rt2  +  3  rt2    .2          3   .2  rt2     - 

+  55  (0,.  +  rl      +  -  09      +30,0.        ri    L      +  -z  0«   0  .  , 

r—  0i  8      Itn        8      2m               ob    lm     4      lm  2m        2      2m  ob) 

kl 


I. ,    -  6     k3     0.0.      0«       *  1  +  ~    r^  (02,    +  §02     +  f  02     ) 
11  —       ob      lm     2m  3     k.        ob        4   lm       4     2m 


(155) 


hi  -  ■  I     ^     «L     *2m   [l  +   10     j2     (0>fa  +  i  »l  ♦  I     0^  )  (156) 

N-  3 


13  N2        ob      lm     2m 


(157) 


k„      ^4 

lm  "2m 


*14  *     I        5      •*-  • 


■30 


".- 


4        N?        2m  k        v   ob  2      lm 


(158) 
(159) 


31  N?  ob        lm        2rn 


(160) 


t  1_1       rt2      .3 

L32   "     4N2        Plm  P2m 


161) 


lk5 


■50  16     N, 


2m 


162) 


'3      ,.2        +  i  0?     +     I  02m  > 


00         N.        ob  k.  ob 


(163) 


+  S     (0\     +I1-     f*       +     1*0*       +502h0?     +^  02,    02  +50*  02b ) 
k1       v   ob  8  lm  8        2m  ob    lm        2  lm  2m        2m  ob' 
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L01 


u .., 


1   +  3 


ih  +i  0L  +  10L> 


(164) 


c   5    ,-4  1    .4  3.4  3    s2   .2  3  „,2  J.  0.2   .2     J 

+  5r—  (0  .    +  q   0.      +  o  0,     +  ^0,0.      +  7  0.   0,     +  30o  0  .    ) 
k  ob        8      1m       8     2m       2     ob   lm       4      lm  2m  2m  ob 


3     k3     rt       M2      1 ;  .  10     k5      eik2  U2         3  .2  . 

X02   s     "2     IT     0ob  0lm    5   *   +     T     k"      (0ob     +  2*1.  +  2  02m> 
1  -  3 


(165) 


r03-      '\    ij     0L     il   +   l0^     <0ob+I0L     +I     0L     >j 


(166) 


"04 


L05 


20 


5  S 

8  N. 
1 

ob 

lm 

k 
1   K5 

16   H 

lm 

3  ^3 
2   N. 

0  , 

0  D 

09 
2m 

10 


k. 


3     „2 


(167) 


(168) 


+  i    0L  40L  >  I       <169) 


2  -'« 


21 


!     if     0lm  0L  L1   +   10     S      (0ob  +  I  0lm     +  1     *L  > 


(170) 


11      k5    «      rt2    <*- 
22  2       -2-     "'ob  *lm  ^2nt 

Nl 


(171) 


-         .       I       S     03       02 
23  4        Nt        lm     *2m 


(172) 


40 


■41 


5 
8 

*5 
Nl 

ob    2m 

5 
8 

k5 

Nl 

0im  0L 

(173) 


(174) 
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